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ABSTRACT 

We present a new method of finding cosmic voids using tomographic maps of Lya 
forest flux. We identify cosmological voids with radii of 2 - 12/i _1 Mpc in a large 
N-body simulation at 2 = 2.5, and characterize the signal of the high-redshift voids 
in density and Lya forest flux. The void properties are similar to what has been 
found at lower redshifts, but they are smaller and have steeper radial density profiles. 
Similarly to what has been found for low-redshift voids, the radial velocity profiles have 
little scatter and agree very well with the linear theory prediction. We run the same 
void finder on an ideal Lya flux field and tomographic reconstructions at various 
spatial samplings. We compare the tomographic map void catalogs to the density 
void catalog and find good agreement even with modest-sized voids (r > 6/i _1 Mpc). 
Using our simple void-finding method, the configuration of the ongoing CL AMATO 
survey covering 1 deg 2 would provide a sample of about 100 high-redshift voids. We 
also provide void-finding forecasts for larger area surveys, and discuss how these void 
samples can be used to test modified gravity models, study high-redshift void galaxies, 
and to make an Alcock-Paczynski measurement. To aid future work in this area, we 
provide public access to our simulation products, catalogs, and sample tomographic 
flux maps. 
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1 INTRODUCTION 

The material we see in the Universe around us makes up 
a beaded, filamentary network known as the “cosmic web” 
(Bond, Kofman V Pogosyan 1996) This web appears to be 
the natural outcome of gravitational instability acting upon 
an initially Gaussian random field. The majority of the cos¬ 
mic web, by volume, is made up of large, almost empty re¬ 
gions known as voids which are surrounded by walls, fila¬ 
ments and clusters (see e.g. van de Weygaert V Platen 2011, 
for a review). Within this paradigm, voids are regions which 
are practically devoid of galaxies. They are slightly prolate 
in shape and occur on a wide range of sizes from Mpc to tens 
of Mpc (Vogeley et al. 1994; Ceccarelli et al. 2006; Lavaux 
& Wandelt 2012). 

The study of cosmic voids has received renewed theo¬ 
retical attention recently. Voids are intrinsically interesting 
as a major constituent of the Universe (by volume) and one 
of the most visually striking features in galaxy maps. They 


* email: caseywstark@berkeley.edu 


form an interesting environment for the study of galaxy evo¬ 
lution. They may present an excellent laboratory for study¬ 
ing material which clusters most weakly (e.g. dark energy or 
massive neutrinos), and for testing modified gravity models. 
Future surveys are expected to find large samples of voids at 
a range of redshifts, enhancing the potential of void science. 

The pristine environments of voids present an interest¬ 
ing setting for the study of early galaxy formation. Galaxies 
in low-redshift voids generally have smaller stellar masses, 
appear bluer, have a later morphological type, and have 
higher specific star formation rates than galaxies in aver¬ 
age density environments (van de Weygaert V Platen 2011; 
Beygu et al. 2015), although the latter properties might be 
solely due to their lower stellar mass (Hoyle et al. 2005; 
Kreckel, Joung V Cen 2011). Extending similar studies to 
higher redshifts to see whether similar trends hold is a press¬ 
ing observational challenge. 

Ryden (1995) was the first to discuss using voids as 
probes of cosmology. Park V Lee (2007) anticipated using 
void ellipticity as a cosmological probe and Lee V Park 
(2009); Bos et al. (2012) discussed constraining dark en- 


© 0000 RAS 



2 Stark et al. 


ergy using voids. Lavaux & Wandelt (2012) investigated 
the potential for using stacked voids as a probe of geomet¬ 
rical distortions (the AP test; Alcock & Paczynski 1979). 
Chan, Hamaus & Desjacques (2014) have studied the clus¬ 
tering of voids and Hamaus et al. (2014b) describe constrain¬ 
ing cosmology with void-galaxy cross-correlations. Hellwing, 
Juszkiewicz & van de Weygaert (2010); Li (2011) have inves¬ 
tigated studying the nature of dark matter using the prop¬ 
erties of voids and Li, Zhao & Koyama (2012); Clampitt, 
Cai & Li (2013) have suggested that void properties may 
provide a strong test of some modified gravity theories. 

Observationally, studies of voids date back over three 
decades (Gregory & Thompson 1978; Longair & Einasto 
1978; Kirshner et al. 1981). Recent redshift surveys have 
identified large samples of voids (e.g. 2dF: Ceccarelli et al. 
2006; SDSS: Sutter et al. 2012, 2014a; VIPERS: Micheletti 
et al. 2014) and a measurement of the AP effect from voids 
in the local Universe has recently been reported by Sutter 
et al. (2014b); Hamaus et al. (2014a). Being underdense in 
both galaxies and dark matter, voids act like objects with 
an effectively negative mass, bending light rays away from 
them. This effect has been recently detected at high signifi¬ 
cance by Clampitt V Jain (2014). 

In the absence of large-scale dynamical and environ¬ 
mental influences, voids would become increasingly isotropic 
as they evolve (Icke 1984). However, in modern theories of 
structure formation the frequent encounters with surround¬ 
ing structures and the influence of large-scale tidal fields 
serve to reverse the simple trend expected for isolated voids 
(van de Weygaert V Platen 2011). As matter in the cen¬ 
ter of voids streams outwards faster than matter towards 
the boundary, the interior evolves into an almost uniform 
low density region surrounded by ‘ridges’ marking the void 
edge: often referred to as a “bucket-shaped” density profile 
(see Ceccarelli et al. 2006; Hamaus, Sutter V Wandelt 2014b, 
for recent fits). The density in the center has a characteristic 
value of S ~ —0.8. 

Historically, surveys of voids over large volumes have 
come from large, galaxy redshift surveys. However, finding 
voids in this manner requires a significant investment in 
telescope time due to the necessity of a high spatial sam¬ 
pling of tracer galaxies. For example, the void catalog pre¬ 
sented in Sutter et al. (2012) found voids in the distribu¬ 
tion of SDSS DR7 galaxies. Their ‘bright’ cut found voids 
with radii larger than 7 /i _1 Mpc with galaxies separated by 
8 /i -1 Mpc. To find comparable galaxy separations at z — 0.5, 
1.0, and 2.0 will require obtaining complete galaxy redshift 
samples for apparent limiting magnitudes of I — 22.5, 24.2, 
and 25.7, respectively (assuming galaxy luminosity functions 
from Dahlen et al. 2005 at z < 1, and Reddy et al. 2008 
at z = 2). So while such galaxy number densities are just 
achievable up to z ~ 1 with existing telescopes, it becomes 
increasingly challenging at higher redshifts. 

In light of the aforementioned challenges, it is under¬ 
standable that little attention has been given to studying 
voids at z > 1 (although see D’Aloisio & Furlanetto 2007; 
Viel, Colberg & Kim 2008). However, recently it has been 
noted that given sufficient sight lines, the Lya forest observed 
in a dense grid of faint background galaxies and quasars 
can be used to create three-dimensional maps of large-scale 
structure and that the observational requirements to map 
out cosmological volumes (V ~ 10 6 h~ 3 Mpc 3 ) are within 


reach of existing facilities (Lee et al. 2014b). Indeed, the 
first pilot map on a small field has already been made us¬ 
ing just a few hours of data from Keck telescope (Lee et al. 
2014a). This method is ideally suited to finding extended 
structures at high redshift. In Stark et al. (2014), we in¬ 
vestigated the possibility of finding protoclusters in tomo¬ 
graphic Lya maps. Here, we study the signature of cosmo¬ 
logical voids in the Lya forest. 

The outline of the paper is as follows. In §2, we briefly 
describe the simulations that we use in this paper and the 
method we use to find voids. In §3, we use our catalog of 
high-redshift voids to explore the properties of voids in den¬ 
sity and flux. In §4, we demonstrate how to find voids using 
tomographic flux maps and discuss how well voids found in 
the maps compare to those found in the matter density. In 
§5, we present the prospects of finding high-redshift voids 
with this method in ongoing and future surveys, and dis¬ 
cuss cosmology applications. We present our conclusions in 
§ 6 . 


2 SIMULATIONS AND VOID FINDING 
2.1 N-body simulations 

In order to study the signal of voids in the Lya forest, we 
make use of cosmological iV-body simulations. We require a 
simulation which simultaneously covers a large cosmological 
volume while having a sufficiently small inter-particle spac¬ 
ing to model transmission in the IGM. The requirements 
are sufficiently demanding that we have used a pure iV-body 
simulation, augmented with the fluctuating Gunn-Peterson 
approximation (FGPA; Petitjean, Mueket & Kates 1995; 
Croft et al. 1998; Meiksin & White 2001; Meiksin 2009). This 
same simulation was also used in Lee et al. (2014b) and Stark 
et al. (2014) so we only review the salient features here. We 
are also providing public access to the relevant simulation 
products - see App. B for more details. 

The simulation employed 2560 3 equal mass (8.6 x 
10 7 H~ 1 Mq) particles in a 256 h ~ x Mpc periodic, cubical box. 
This provides sufficient mass resolution to model the large- 
scale features in the IGM at z = 2-3 using the FGPA 
(Meiksin & White 2001; Rorai, Hennawi & White 2013) and 
sufficient volume to find large voids. The assumed cosmology 
was of the flat ACDM family, with Q m ~ 0.31, Ub^ 2 ~ 0.022, 
h = 0.6777, n s = 0.9611, and as = 0.83, in agreement with 
Planck Collaboration et al. (2013). The initial conditions 
were generated using second-order Lagrangian perturbation 
theory at Z[ c — 150, when the rms particle displacement was 
40 per cent of the mean inter-particle spacing. The parti¬ 
cle positions and velocities were evolved using the TreePM 
code of White (2002). Throughout the text, we will use the 
particle positions and velocities from the output at z = 2.5. 
Using the particle positions and velocities at z — 2.5, we 
generated mock Lya forest spectra on a 2560 3 grid with the 
FGPA as described in Stark et al. (2014). In all, we gener¬ 
ated 2560 3 grids with the matter density and Lya forest flux 
in real- and redshift-space. For many purposes in this work, 
we did not need the high resolution provided by the 2560 3 
grids and found it much easier to work with smaller grids. 
For this reason, we also downsampled the fields to 256 3 by 
simply averaging neighboring grid points. In the remainder 
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of the paper, when we refer to flux, we mean the Lya forest 
transmitted flux fraction perturbation 

$f = F/(F) - 1 . (1) 

2.2 Void finding 

There are a variety of methods and tools used to find voids 
in large-scale structure (Kauffmann & Fairall 1991; Platen, 
van de Weygaert & Jones 2007; Neyrinck 2008; Sutter et ah 
2015). We use a simple spherical underdensity method on 
the low resolution, gridded densities to construct our z = 2.5 
void catalog. This technique is similar to spherical overden¬ 
sity (SO) halo finding, but instead applied to underdensities. 
We identify voids by taking the 256 3 density grid and select¬ 
ing points under some threshold value, then growing spheres 
around the points until the average density enclosed reaches 
a target value. We handle overlapping voids by only saving 
the void with the largest radius, and we also discard any 
remaining voids with a radius r < 2 h -1 Mpc. 

The free parameters in this void-finding method are the 
threshold value pthresh and enclosed average target value 
Penc • For the density field, there is a well-motivated thresh¬ 
old value of pthresh = 0.2 p. This threshold value is a canon¬ 
ical density for a void core, corresponding to the central 
density at shell-crossing for spherical void models (van de 
Weygaert & Platen 2011). This critical density value is a 
standard choice in other void-finding codes (e.g. Neyrinck 
2008). The choice of the average target value, however, is 
somewhat arbitrary. We first tested an average target value 
of penc = 0.2 p, but found that it produced voids that were 
far too small - the spheres never reached the apparent ‘edge’ 
surrounding the low-density core. We experimented with 
several other average target values and found that a value of 
Penc = 0.4 p resulted in good agreement between the sphere 
sizes and the apparent void edges. 

In principle, there is nothing special about the specific 
threshold and average density values we chose, and these 
parameters should depend on the redshift. That is, as voids 
continue to evacuate, the core and average densities of the 
void will decrease. In practice, we found the final void cata¬ 
log is not very sensitive to these settings, although the void 
radii clearly scale with the average target density setting. 
Since most large voids have central densities < 0.2 p already, 
the exact value of the threshold mostly makes a difference 
in terms of how many points we must search over, and less 
of a difference in the void centers. We did find that a very 
small threshold (say < 0.1 p for this redshift) forces voids 
to grow from positions that often look off-center by eye. 
This is due to the hierarchical nature of voids, in that the 
lowest-density point in a large void is typically the center 
of a smaller subvoid, sometimes referred to as the void-in¬ 
void scenario (c.f. fig. 6 of Neyrinck 2008 and Sheth & van 
de Weygaert 2004). Using the SO parameter values of the 
threshold pthresh = 0.2 p and average target p enc = 0.4 p, we 
found 16,167 voids which cover 15 per cent of the simulation 
volume. 

Fig. 1 shows a slice through our simulation, centered on 
one of the largest voids in our catalog with r — 9.7 /i _1 Mpc. 
The slice is 40/i _1 Mpc across and 6/i _1 Mpc projected into 
the page. The four panels show the void in real-space 


density, redshift-space density, Lya forest flux, and a to¬ 
mographic flux map constructed from a mock survey. In 
each image, we marked the void center and radius with a 
black dot and dashed line, respectively. In the first panel, 
we also overplotted the positions of halos found in the 
same slice. Green dots mark the positions of halos with 
mass M ^ 10 12 (roughly an L* halo at this red¬ 
shift), while yellow dots mark the positions of halos with 
3 x 10 11 h~ 1 M & ^ M < 10 12 Based on simple abun¬ 

dance matching (see Fig. 2), these halos should host galaxies 
with apparent magnitudes 1Z < 24.7 and 24.7 < 1Z < 25.6, 
just bright enough for redshift determination with existing 
facilities. The second panel shows that in redshift space the 
void has a larger density contrast and extent in the line- 
of-sight direction due to the outflow of matter from the 
void. Such a large structure is easily visible in the redshift- 
space density and flux. Although the tomographic flux map 
is a blurred version of the true flux, the void structure is 
so large that it can still easily be picked out by eye. For 
reference, the tomographic map is one of the realizations 
from Stark et al. (2014) with an average sightline spacing 
of d± = 2.5h -1 Mpc, similar to the ongoing survey of Lee 
et al. (2014a). At the same time, the void is not captured 
by the galaxy positions even if we assume a complete galaxy 
sample. The relative sparsity of such halos highlights the 
difficulty in finding voids, even large ones, at high redshift 
using galaxies as tracers. 

We compared our void catalog with that produced by 
a watershed void finder, similar to the ZOBOV code. The 
watershed method finds the set of connected elements all 
under some threshold. In ZOBOV the elements are the 
Voronoi cells in the tessellation of the dark matter particle 
positions (where the density is estimated from the volume 
of the Voronoi cell), but in this case, we use the density 
values on the 2560 3 grid for simplicity. The watershed al¬ 
gorithm on a uniform grid is straightforward. We find the 
set of points under the given threshold, and keep a list of 
the under-threshold points that have not been assigned to a 
specific watershed. Starting with the minimum value point, 
we search grid neighbors to see if they are also under the 
threshold and add them to the current watershed if so. The 
search stops when there are no remaining neighbors under 
the threshold. These points are then removed from the list 
of unassigned points and we move on to the next watershed. 
After we assign all points under the threshold, we discard 
watersheds with an effective radius r e ff = {3V s h.ed/^^ 
less than 2 /i _1 Mpc, as we did with the spherical underden¬ 
sity voids. Using this method with the same threshold of 
p < 0.2 p, we found 6,364 voids, covering 5 per cent of the 
simulation volume. 

The sets of large voids in the spherical underdensity (or 
SO) catalog and the watershed catalog agree very well. We 
visually inspected the 100 largest voids in the SO catalog, 
and found matches in the watershed catalog. In most cases 
the watershed void effective radius was slightly smaller (by 
l-2/i -1 Mpc), which explains the total count and volume 
difference, and the watershed voids typically have complex 
morphologies. The watershed voids often have an ellipsoidal 
core, with fingers stretching out into smaller low-density re¬ 
gions. We compared the SO void centers to the watershed 
void value-weighted centroids x s h e d = ^i^ip^ 1 / JV p” 1 , 
where the sums are over all the points in the shed, and we 
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Figure 1. A slice through our simulation and a mock reconstructed flux map centered on a large void, with a radius of 9.7h -1 Mpc. 
The slice is 40h -1 Mpc across and 6h -1 Mpc into the page. The vertical axis, xp is the redshift direction, along the line of sight, and 
the horizontal axis, x±, is one of the directions transverse to the line of sight (the other transverse direction is into the page). (Left) 
The matter density, in real space. The void center and radius are indicated with a black dot and dashed circle. Positions of halos with 
M ^ 10 12 h~ 1 M@ are plotted as green dots and positions of halos with M ^ 3 x 10 11 are plotted as yellow dots. These halos 

would host galaxies at the spectroscopic limits for existing instrumentation. (Middle left) The matter density, in redshift space. (Middle 
right) The ‘true’ flux field over the same volume, in redshift space. (Right) A reconstructed flux field from a mock survey (also in redshift 
space) with average sightline spacing, d± = 2.5h -1 Mpc (see text). 


weight by the inverse of the density so that the centering 
is driven by lower-density points. Unfortunately, the non- 
spherical geometries of the watershed voids tend to drive the 
centroid away from the center found with the SO method 
and the centers in the two catalogs tend to disagree by sev¬ 
eral Mpc (see Appendix A for more discussion and images of 
the watershed voids). It is reassuring that these two methods 
for finding voids in the density field qualitatively agree well, 
but we decided to use the SO void catalog for the remainder 
of this work due to its simplicity. Overall, the centers and 
simple shapes of the SO voids provide cleaner radial pro¬ 
files and should be easier to find in the tomographic maps 
later. We were also concerned that the non-trivial noise we 
expect in the maps from tomographic reconstruction might 
artificially combine or split watershed regions, whereas the 
spherical average in the SO method will be less affected by 
such noise. Since our tomographic maps come with a noise 
estimate, one could imagine a more sophisticated algorithm 
(e.g. a matched filter or likelihood-based method) for finding 
voids could be implemented. We leave such investigations to 
future work. 


3 VOIDS AT z = 2.5 

In Fig. 3, we show the cumulative number density of voids 
as a function of void radius. We plot the distribution of our 
z = 2.5 voids in black, and show the distribution of low- 
redshift voids with a green dashed line (z ~ 0.5), computed 
from eq. 21 of Lavaux & Wandelt (2012). As expected, there 
are many more small voids and voids are generally smaller 
at z = 2.5 than at z ~ 0 (c.f. fig. 1 of Ceccarelli et al. 2006 
or fig. 7 of Lavaux &; Wandelt 2012). While voids with radii 
of 7h -1 Mpc are common for low-redshift studies, we have 
only 126 voids with r ^ 7 /i _1 Mpc, which cover two per 
cent of the simulation volume. We note, however, that it is 
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Figure 2. The cumulative, halo mass function in the simula¬ 
tion at z = 2.5, as a function of (FoF) halo mass. The dotted 
lines indicate the abundances of galaxies brighter than the listed 
apparent TUmagnitude limits, derived from the luminosity func¬ 
tion of Reddy et al. (2008). Since at this redshift and these masses 
satellites make up a small fraction of galaxies by number, this plot 
allows us to approximately equate our mass-limited halo catalogs 
into flux-limited galaxy catalogs. Note we have used volumes and 
masses without factors of h in this figure to match the scalings 
adopted in Reddy et al. (2008). 

difficult to compare void sizes across works using different 
void-finding methods and working at different redshifts. For 
instance, we could increase the number of r ^ 7 h x Mpc 
voids by simply increasing the average target value in our SO 
void finder. For the most part, this does not change which 
large voids are identified, but does shift centers and increase 
the cumulative number density at a particular value. See 
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Figure 3. The cumulative, comoving number density of voids 
with radii greater than r vs. r. In black, we show the distribution 
of our SO density voids at z = 2.5. The green dashed line is 
the distribution of the voids found in Lavaux & Wandelt (2012), 
which are generally larger, as expected. 

Colberg et al. (2008) for more detail about the difficulties of 
defining voids and differing results from various void finders. 

Fig. 4 shows the radial profile of voids in density, radial 
peculiar velocity, and real-space flux, stacked by radius. The 
gray lines show the profile for individual voids, while the 
thick black lines show the mean. We chose to stack the voids 
in our catalog with radii in the range 5 ^ r < 6/i _1 Mpc, 
of which there are 511. We note that the central density is 
about 0.2 p in the dark matter distribution, as a result of our 
choice of void finder. The average profile then rises almost 
continuously to the void edge, though individual voids show 
substructure within them (also visible in Fig. 1). The slope 
of our profile contrasts with the profile of more evolved voids 
at lower redshifts, which exhibit a ‘bucket’ profile 1 . We used 
the z — 0 simulation output to create a low-redshift SO void 
catalog (with adjusted threshold and average target values) 
and found that these voids do exhibit such a ‘bucket’ profile. 
The voids we are studying at z — 2.5, however, have not 
yet evolved to such a state and are still in the process of 
evacuating. The real-space flux profiles illustrate just how 
well the flux profile mirrors the density profile. The center 
of the voids have F ~ 1, which translates to Sf — 0.25 
for our setting of ( F) = 0.8. In both the density and flux 
value, the stack profile almost reaches the mean value by 
r — 10 h~ x Mpc. 

The middle panel of Fig. 4 shows the radial velocity pro¬ 
files of the voids. The profiles are linearly increasing up to 
the void radius, where they peak around 120kms _1 before 
turning over. The radial velocity profiles also have a fairly 
small scatter - at the average radius of 5.4/i _1 Mpc, the 
mean velocity is 119.7kms -1 with a standard deviation of 
5.4kms _1 or about 5 per cent. Within the context of linear 
theory v(f) oc VV -2 5(r). If we assume a spherical mass dis¬ 
tribution, this can be solved to yield v r oc S(< r)/r 2 , where 

1 For example Ceccarelli et al. (2006) propose p(r)/p = Aq + 
A 3 ( r/Ry ) 3 while Hamaus, Sutter & Wandelt (2014a) proposes 
a 4-parameter model with a similar shape. 




Figure 4. Radial profiles centered on voids with 5 ^ r < 
6/i -1 Mpc. Gray lines are individual profiles, and the thick black 
line shows the mean. (Top) The matter density profile. (Middle) 
The radial peculiar velocity profile. The dashed red line shows 
the linear theory prediction using the density stack profile above. 
(Bottom) The real-space flux profile. 

S(< r) is the “overdensity enclosed within r” in analogy with 
the Newtonian gravitational acceleration due to a spherical 
mass distribution. The radial velocity profile around spheri¬ 
cal, or averaged, voids then becomes (Peebles 1993; Hamaus, 
Sutter & Wandelt 2014a) 2 

v(r) = -afH± £ - l) x 2 dx (2) 

where / ~ ^m 55 (^) is the growth factor, which at z = 2.5 
is close to 1. This form was shown in fig. 1 of Hamaus, Sut¬ 
ter & Wandelt (2014a) to fit the velocity profile of stacked 
voids in N-body simulations well at z ^ 0. The dashed red 
line in Fig. 4 shows this linear theory approximation, which 
we see compares favorably to the profile measured in our 
z = 2.5 voids (within 10 per cent over the range plotted). 

2 Note, we have an additional factor of a in this expression com¬ 
pared to Equation 2.2 of Hamaus, Sutter & Wandelt (2014a), 
perhaps due to a difference in proper vs. comoving quantities. We 
always use comoving scales and densities, and peculiar velocities. 
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Figure 5. Contour plots of a stack of the voids with 5 ^ r < 
6h -1 Mpc, showing the impact of redshift-space distortions. The 
four panels show the density and flux in real and redshift space, 
binned in the distances parallel (xy) and perpendicular (x±) to 
the line of sight. In the redshift-space fields, the radial profiles are 
clearly extended in the line-of-sight direction. 


It is somewhat surprising that the linear theory prediction 
matches our simulated radial velocity profile result down to 
Mpc scales and for \S\ = 0.8. The fact that this prediction 
also matched void radial velocity profiles at z = 0, with 
voids from a different finder method is impressive (Hamaus, 
Sutter & Wandelt 2014a). 

Fig. 5 shows the two dimensional profiles (in mass and 
flux) of stacks of voids with radii 5 ^ r < 6 A -1 Mpc in 
both real and redshift space. Apart from some noise near the 
line-of-sight axes, the contours in Fig. 5 are isotropic in the 
real-space panels but show extended, anisotropic profiles in 
redshift-space. This is an indication of the effect of peculiar 
velocities, which appear visually to be larger in our case than 
at lower redshifts when voids are traced by galaxies. We note 
that the profiles are much better measured at small radii 
where there is less scatter. Beyond the radius of the stack r ~ 
5.5 A -1 Mpc, the scatter in the individual profiles increases 
significantly. We believe this is the source of the extended 
orange contour in the bottom right panel, for instance. 

Since we expect the stacked voids to be isotropic in real 
space, by symmetry, any observed anisotropy offers an op¬ 
portunity to study such peculiar velocities. This could be 
particularly interesting for constraining models with mod¬ 
ified gravity. For example, Clampitt, Cai & Li (2013) find 
that, driven by the outward-pointing fifth force, individual 
voids in chameleon models expand faster and grow larger 
than in a ACDM universe. Such effects would modify the 
profile of the stacked voids in a potentially observable man- 


Table 1. Void catalogs 


Field 

SO Thresh. 

SO Avg. 

Count 

Vol. Frac. 

P 

0.20 p 

0.4 p 

16,167 

0.152 

Pred 

0.15 p 

0.3 p 

16,338 

0.151 

8 p 

0.224 

0.167 

16,296 

0.150 

hires map 

0.224 

0.167 

16,586 

0.203 

midres map 

0.224 

0.167 

8,724 

0.181 

lores map 

0.224 

0.167 

5,565 

0.153 


The spherical under/overdensity void catalogs used for compari¬ 
son, found in our V ~ 1.7 X 10 7 h -3 Mpc 3 simulation. We use the 
original p catalog as our ‘truth’ and varied the SO parameters 
for the p rec j and 8 p catalogs to qualitatively match. The hires 
map uses a d = 2.5 h~ 1 Mpc sightline spacing mock survey, while 
the midres map comes from a d = 4.0h —1 Mpc configuration, and 
the lores map from a d = 6.0h -1 Mpc configuration. The noise 
and smoothing inherent in the tomographic reconstruction pro¬ 
cess create larger differences in the catalog properties. 


ner, allowing observations of voids in the Ly<a forest to test 
such models. Based on the radial velocity profiles shown 
in Fig. 4 and the measured standard deviation, one would 
need only about 20 voids with a radial velocity measurement 
to reach one per cent standard error (assuming Poisson er¬ 
rors). With accurate enough radial velocity measurements 
from void anisotropies, it should be possible to detect devi¬ 
ations at the 10 per cent level with relatively small samples. 
Conversely, the larger impact of redshift-space distortions in 
the Lyce flux field means they must be modeled in order to 
make a measurement of the Alcock-Paczynski effect (Alcock 
& Paczynski 1979) from stacked voids in the flux field (see 
Section 5.2). 


4 FINDING VOIDS IN FLUX 

Underdense regions show up as high transmission regions in 
the 3D Ly<a forest flux for z = 2 - 3, as shown in Fig. 1 and 
4. This is not necessarily the case at lower redshifts, since 
the characteristic density probed by the forest increases with 
time (Becker et al. 2011; Lukic et al. 2015). At lower red- 
shifts, it is difficult to see differences in transmission passing 
through an underdense region vs. a moderately overdense 
region, since it takes a significant overdensity to create an 
observable absorption feature. Fortunately, there is a large 
overlap between the redshift range of the forest accessible 
from the ground and the redshift range where mean den¬ 
sity structures scatter an observable fraction of the light. 
Given this, finding voids in flux at z = 2 - 3 is a matter 
of finding coherent high-transmission regions. In Stark et al. 
(2014), we outlined a simple method using tomographic flux 
maps to find protoclusters (coherent low-transmission re¬ 
gions) at these redshifts. We now adapt these methods to 
find coherent high-transmission regions, corresponding to 
high-redshift voids. 

The void catalogs used in this section are listed in Ta¬ 
ble 1, including voids found in the redshift-space density, 
the flux, and three tomographic flux maps. In each case, we 
modified the SO threshold and average target parameters to 
create a void catalog with roughly the same void count and 


© 0000 RAS, MNRAS 000 , 1-14 















Voids in Lya Forest Tomographic Maps 7 


radius distribution as the real-space density catalog. For ref¬ 
erence, we use three of the tomographic flux maps created in 
Stark et al. (2014). These maps were constructed by mock¬ 
ing up a realistic survey covering the simulation volume with 
signal-to-noise distributions similar to the pilot observations 
of Lee et al. (2014a), and several other settings of the aver¬ 
age sightline spacing and minimum spectral signal-to-noise 
ratio. The tomographic maps were then generated by run¬ 
ning our reconstruction code on the mock spectra. The three 
maps we use are generated from mock surveys with aver¬ 
age sightline spacings of d± = 2.5, 4, and 6h _1 Mpc. The 
smallest sightline separation configuration is similar to the 
ongoing CLAMATO survey, and the larger separation con¬ 
figurations are similar to what we expect from large-area 
surveys on 8 - 10 m telescopes like the Subaru Prime Focus 
Spectrograph (PFS; Takada et al. 2014) or the Maunakea 
Spectroscopic Explorer (Simons et al. 2014). We refer to 
these tomographic maps as the hires, midres, and lores flux 
maps. We discuss the characteristics of the individual cata¬ 
logs in the following subsections and focus on our method to 
compare catalogs for now. Comparing across void catalogs, 
we want to confirm that there are nearby pairs of voids with 
similar radii. We expect to find the same set of voids in den¬ 
sity and flux, as we have demonstrated how well-matched 
the fields are in previous sections. This is mostly a matter 
of determining the best SO parameters for the flux. The 
tomographic flux maps, however, are contaminated by and 
spectral noise in individual mock spectra and shot noise due 
to sparse sampling of the field, and this noise will certainly 
affect our capability of finding voids. 

We use two metrics to compare the catalogs of voids 
found in different fields. The first metric is essentially the 
sum of the difference in the center positions and radii, which 
we call the match error. If we are comparing voids in catalog 
A with voids in catalog B, for each A-B pair, we compute 
the match error 

_ V Jta ~ r B ) 2 + |XA - X B | 2 /3 2 

VA 

where r is the radius and x is the center position. We chose 
this form of the error for several reasons. First, this form of 
the error also allows for trading off differences in radii and 
centers. We want to consider the differences relative to the 
size of the void, which will allow for larger center and radii 
differences for larger voids. Note that this form of the error 
assumes the radius of void A is the reference. Finally, we 
compare 1/3 of the center difference to the radius difference 
just due to the dimensionality (and empirically we found 
that the mean center difference is about 3 times the mean 
radius difference). Later in this section, we will show that a 
match error e < 0.3 qualifies as a good match for a void be¬ 
tween two catalogs, and we will use this cut to count which 
voids are ‘matched’. The second metric we use is the total 
volume overlap between the voids in both catalogs. Clearly, 
this metric is less useful for telling if a catalog A void is well- 
matched by a single catalog B void. However, it is a useful 
measure of how well-matched the catalogs are overall and 
does not depend on a specific form of the error nor a spe¬ 
cific value to cut at. It is also useful in cases where noise in 
the tomographic maps artificially combines or splits voids - 
although the centers and radii might not match across cata¬ 
logs, there will still be a sizable volume overlap. We can use 


these metrics to get a sense of void completeness and purity 
of each of the flux catalogs with respect to the density void 
catalog. The number of ‘matched’ density voids compared to 
the total number of density voids (the match fraction) is a 
measure of completeness. We also measure completeness by 
comparing the overlapping volume between two catalogs to 
the total volume in density voids (the overlap fraction). The 
purity of the flux catalogs can be measured by matching in 
the other direction (the fraction of matched flux voids) and 
by comparing the overlap to the total volume in flux voids. 

Using these metrics, we first found that redshift-space 
distortions can create large differences in the centers and, 
to a lesser extent, the radii of the voids. In order to more 
easily compare voids found in density and flux, we created 
a void catalog using the redshift-space density. Voids found 
in redshift-space density matched those found in real-space 
density best (in terms of detecting the same voids with sim¬ 
ilar radii) when we used a threshold of 0.15 p and an average 
target of 0.3 p. These densities are lower than the real-space 
values since outflows from voids drive densities lower. 

Before applying these metrics to the void catalogs de¬ 
rived from the various flux maps, we also compared the 
redshift-space density void catalog to random void catalogs, 
mainly to get a sense of the worst-case performance. We cre¬ 
ated ten catalogs of 16,338 voids (the same number as the 
redshift-space density catalog), with centers uniformly dis¬ 
tributed in the simulation domain, and with radii randomly 
drawn from the same distribution as that in the redshift- 
space density catalog. We compared each random void cat¬ 
alog against the density catalog, computing the fraction of 
density voids with a match error e < 0.3 and the fraction of 
the total volume overlap to the total volume in density voids. 
Overall, 2.7 per cent of the density voids were matched by 
voids in the random catalogs on average. It is reassuring to 
see that a small fraction of the density voids are matched 
by random voids which tells us that our cut of e < 0.3 is 
stringent enough. We also noticed that for the largest voids 
(r ^ 8h _1 Mpc), the average match fraction drops to 1.3 
per cent. This is due to the fact that both the density and 
random catalogs contain just a few very large voids and it 
is even less likely that they will overlap enough to meet the 
match error cut. The average volume overlap fraction be¬ 
tween density voids and voids in the random catalogs was 
15 per cent, and did not change with the radius considered. 
This is not surprising since the voids cover roughly 15 per 
cent of the total volume, so random points will overlap about 
that often. 

4.1 Ideal flux 

We ran our SO void finder on the ‘true’ flux grid (256 3 ) 
using a threshold of Sf ^ 0.224 and an average target of 
Sf — 0.167. We experimented with several values of the 
threshold and average target fluxes and found that these 
values resulted in a number of voids and radius distribution 
similar to the catalog of voids found in density. The mapping 
from flux to density evolves quickly with redshift, so these 
SO parameters would have to be adjusted for other redshifts 
and UV background prescriptions. 

The flux void catalog matches the redshift-space density 
catalog very well. For all voids (r ^ 2 h -1 Mpc), 84 per cent 
of the density voids and the flux voids are matched. The 
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Table 2. Catalog comparison for r ^ 6 h 1 Mpc voids 



Density 

Flux 

Hires map 

Midres map 

Lores map 

Random 

Density 

Flux 

Hires map 
Midres map 
Lores map 
Random 

0.988 / 0.933 
0.581 / 0.567 
0.284 / 0.409 
0.186 / 0.349 
0.0193 / 0.151 

0.994 / 0.937 

0.576 / 0.567 
0.282 / 0.408 
0.184 / 0.347 
0.0191 / 0.150 

0.660 / 0.631 
0.683 / 0.637 

0.258 / 0.425 
0.191 / 0.377 
0.0368 / 0.204 

0.478 / 0.568 
0.471 / 0.569 
0.356 / 0.484 

0.192 / 0.359 
0.0480 / 0.182 

0.269 / 0.517 
0.292 / 0.514 
0.238 / 0.430 
0.171 / 0.344 

0.0398 / 0.149 

0.0194 / 0.152 
0.0188 / 0.153 
0.0261 / 0.153 
0.0243 / 0.152 
0.0204 / 0.153 


The catalog void match fraction and volume overlap fractions (separated by a slash in each cell), for r ^ 6h -1 Mpc voids. For each row, 
we compute the fraction of voids with a match error e < 0.3 and the fraction of the total volume overlapped by voids in the catalog of 
that column. 


volume overlap fraction is also very high, at 86 per cent of 
the density void volume and also 86 per cent of the flux void 
volume. For larger voids (r ^ 6/i _1 Mpc), the catalogs are 
even better matched. In this case, 99.4 per cent of the density 
voids and 98.8 per cent of the flux voids are matched, and 
the volume overlap fraction is 93.7 per cent of the density 
voids and 93.3 per cent of the flux voids. For reference, the 
density catalog contains 335 of these large voids, while the 
flux catalog contains 325. 

4.2 Tomographic flux maps 

We constructed the map void catalogs by running the SO 
void finder on the maps with the same SO parameters we 
used for the ideal flux field. We tried several other SO pa¬ 
rameter settings on the maps, but found that the default 
value catalogs performed best in comparison to the den¬ 
sity and flux catalogs. Small changes to the SO parameters 
resulted in slightly better performance, but changes larger 
than about A Sf = 0.01, resulted in similar or worse per¬ 
formance, so we did not bother optimizing these parame¬ 
ter choices further. Unfortunately, our simple void-finding 
method does not consider noise in the map which can con¬ 
taminate the set of thresholded points and the spherical 
averages. The noise in the map acts to scatter points be¬ 
low or above the SO threshold, creating false negatives and 
positives respectively. In the same way, the noise can af¬ 
fect the spherical averages used in the SO finder, resulting 
in inaccurate radii. However, this should be less of an is¬ 
sue for coherent structures spanning several map resolution 
scales, which is apparent in our results for small vs. large 
voids. The effects of the noise are apparent in the radius 
distribution of the map catalogs. In the hires map cata¬ 
log, the number of very small voids (r < 2.5/i _1 Mpc) is 
3,902, about half of the number found in the density cata¬ 
log (6,157). This is likely due to shot noise where sightlines 
did not sample these smaller structures well enough. The 
number of medium voids (3/i _1 Mpc < r < 6/i _1 Mpc) is 
about double that in the density catalog, and the number 
of large voids (r ^ 7 h~ 1 Mpc) about the same (147 vs. 121). 
This explains why there is a similar total number of voids in 
the density and hires map catalogs, but more total volume 
in the map catalog (see Table 1). The radius distributions 
of the midres and lores map catalogs are more distorted by 
the noise. The midres map catalog contains about a half 
the number of small voids (r < 5/i _1 Mpc) compared to 



2-4 4-6 6-8 >8 

r bin ( h~ l Mpc) 


Figure 6. A box and whisker plot of the distributions of match 
errors vs. radius. The match error is computed between the 
redshift-space density voids and the voids in the flux catalog 
(gray), the hires map catalog (blue), and in one of the random 
catalogs (red). The line in the middle of the box is the median, the 
box extends from the 25th to the 75th percentile, the whiskers ex¬ 
tend down to the 5th and up to the 95th percentiles, and points 
outside of this range are plotted individually. The match error 
threshold of e = 0.3 is marked with a dashed line. The flux cata¬ 
log matches the density catalog exceptionally well. The hires map 
catalog is essentially random for small voids, but performs much 
better for large voids. The random catalog match error is fairly 
flat across radius bins. 


the density catalog and the lores map catalog contain about 
a quarter. There are approximately double the number of 
large voids (r ^ 7 h _1 Mpc) in both map catalogs compared 
to the density catalog. 

In Fig. 6, we plot distributions of the match errors be¬ 
tween the voids in the redshift-space density catalog and 
voids in the flux, hires map, and one of the random catalogs 
vs. radius. The line in the middle of the box shows the me¬ 
dian, the box extends from the 25th to the 75th percentiles, 
the whiskers extend from the 5th to the 95th percentile, and 
samples outside this are plotted individually. The match er¬ 
ror cut value of e = 0.3 is marked with a dashed line. The 
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match errors against the random catalog tend to fall around 
e = 0.6, and there are few points under our cut of e = 0.3, 
again showing that this is a safe choice. It is also reassuring 
that the random errors are relatively flat over radius bins be¬ 
cause we defined the match error relative to the original void 
radius. The gray distributions show just how well-matched 
the flux and density catalogs are and that the flux match er¬ 
rors overlap very little with the random errors. Overall, the 
hires map catalog misses a significant fraction of the small 
voids in the density catalog, but performs well for larger 
voids. For all voids (r ^ 2/i _1 Mpc), the hires map catalog 
matches only 16.3 per cent of the density voids and the vol¬ 
ume overlap fraction is 49.4 per cent. This can also be seen 
in the smallest radius bin in Fig. 6, where the hires map 
is just a bit lower than the random distribution. However, 
considering larger voids (r ^ 6/i _1 Mpc), the hires map cat¬ 
alog performs much better matching against 66 per cent of 
the density voids and overlapping with 63 per cent of the 
volume. In Fig. 6, there is a clear trend that the hires map 
match errors decrease with radius, separating from the ran¬ 
dom distribution. The match fractions in the other direction 
(purity of the hires map voids) are similar at 17.8 per cent 
for all voids and 58.0 per cent for large voids. The lower 
match fractions for large voids in this case is driven by the 
hires map catalog having more large voids. 

The midres map void catalog performs worse for all 
voids, but still matches a considerable fraction of the den¬ 
sity voids. Overall, the midres map catalog matches only 6.2 
per cent of the density voids, although it still overlaps with 
40 per cent of the density void volume. If we consider larger 
voids (r ^ 6h -1 Mpc), the midres map catalog matches 48 
per cent of the density voids, and matches 60 per cent of 
even larger voids (r ^ 8/i -1 Mpc). The lores map catalog 
performs worse than this, but is still useful for finding large 
voids. The lores map catalog matches only 2.5 per cent of 
the density voids - consistent with the random catalogs - 
although it overlaps with 30 per cent of the volume. This 
indicates that the void finder is still able to find regions 
containing voids from the map, but does not recover an ac¬ 
curate center or radius. If we consider some of the largest 
voids (r ^ 8h _1 Mpc), the match fraction increases to 48.9 
per cent and the volume overlap fraction is 61.9 per cent, 
again confirming that the larger the void, the better the 
maps perform. 

Our results are also summarized in Table 2. In this ta¬ 
ble, we give the match and overlap fractions between the 
redshift-space density, flux, hires map, midres map, lores 
map, and random catalogs for the voids with r > 6 h 1 Mpc. 
We note again that the first five catalogs are single catalogs 
while the random results are averages over the ten random 
catalog realizations. The trends between the catalogs are the 
same as described above: the random catalogs match 2-3 
per cent and overlap about 15 per cent, and the correspon¬ 
dence between the density and flux catalogs is very high. 
Comparing the density (or flux) catalogs to the map cata¬ 
logs, the match fraction drops to 60 - 70 per cent for the 
hires case, to 40 - 50 per cent for the midres case, and down 
to 20 - 30 per cent for the lores case. However, the volume 
overlap fraction remains relatively high for all of the maps 
indicating that the poor matching is more the fault of our 
simple void finding method than the maps truly missing the 
voids. 


Overall, the maps perform decently matching voids with 
radii larger than the map resolution, but it is surprising that 
the maps still do not perform better for the largest voids. 
We visually inspected many of the large voids to see why 
the flux map void catalogs sometimes miss these large voids. 
We show two example slices of large voids in the flux, hires 
map, and lores map in Fig. 7. The top row shows a suc¬ 
cessful void match in both maps, while bottom row shows a 
failure case. In the flux panels, we overplot the original void 
with a black dot and circle. In the map panels, we overplot 
the best match void and annotate the match error. In the 
top panels, the matching hires map void is a bit smaller 
and offset just a bit to the bottom left. In the lores map, 
a noise feature around (x,y) = (5, —5) h _1 Mpc pushes the 
matching void center up farther, but with a similar radius, 
resulting in a sufficiently small match error. In the bottom 
hires map panel, there is significant noise around the cen¬ 
ter of the original void which pushes the center of the void 
up and restricts the growth of the void radius to a much 
smaller size. In the bottom lores map panel, the filamentary 
structure around (x,y) = (5, — 5) h _1 Mpc is missing, which 
allows the void radius to grow much larger and results in 
a poor match. Interestingly, in the bottom row, the over¬ 
all structure of the lores map matches the structure of the 
ideal flux better than the hires map by eye. However, we 
have not considered estimate of the noise in the map. Using 
the tomographic reconstruction method outlined in Stark 
et al. (2014), it is possible to compute the covariance of the 
map or to run Monte Carlo error estimates. Considering the 
amount of noise apparent in Fig. 7, much could be gained 
by incorporating a noise estimate into a void finding pro¬ 
cedure. We believe future work can make significant gains 
in void finding performance by considering the structure of 
voids beyond simple spheres and taking the map noise into 
account. 


5 DISCUSSION 

5.1 Survey prospects 

For the cosmology of our simulation, the comoving ra¬ 
dial distance to z = 2.5 is 4050/i _1 Mpc, thus one degree 
subtends 70/i _1 Mpc. Assuming a 250/i _1 Mpc depth (e.g. 

2.2 < z < 2.5), each square degree of survey area translates 
into a volume of 1.2 x 10 6 /i -3 Mpc -3 . Given the number 
densities in Fig. 3, we see that surveys like CLAMATO with 
V ~ 10 6 /i _3 Mpc 3 would encompass about 150 voids larger 
than 5/i -1 Mpc in radius. If we assume a conservative void 
finding efficiency of 60 per cent, our simple method would 
recover ~ 90 voids. This would be the first detection of a 
significant population of high-redshift voids. Of course, this 
is a lower bound on the efficiency of identifying voids with 
a map of this resolution due to our conservative choice of 
what constitutes a match, and that there is still room for 
improvement in the method. Using the Subaru Prime Focus 
Spectrograph (PFS; Takada et al. 2014), it is possible to 
double the target density, covering a larger redshift range at 
the cost of sightline density. In Lee et al. (2014b), we dis¬ 
cussed piggybacking on the planned galaxy evolution sur¬ 
vey described in Takada et al. (2014). Such a survey would 
provide a map of roughly 16 deg 2 or 8 x 10 4 h _2 Mpc 2 area 
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Figure 7. Slices of two large voids showing the flux, hires map, and lores map (from left to right). A void that is well-matched in the 
map void catalogs is shown on top, and a poor match case is shown on the bottom. The top void has a radius of 10.9h -1 Mpc and the 
bottom void has a radius of 8.3h -1 Mpc. Just as in Fig. 1, the slice is 40h -1 Mpc across and 6h -1 Mpc into the page, although in this 
case, the redshift direction is into the page. The black dot and circle in each panel show the void as found in each catalog. The flux 
panels show the original void (very close to the one found in density), while the map panels show the map voids with the lowest match 
error. We also annotate the match error of the map voids in the top-left corner of the images. 


and 700/i -1 Mpc depth (2.3 < z < 3.2) for a total vol¬ 
ume of 6 x 10 7 /i _3 Mpc 3 , although at a coarser resolution of 
about 5 /i -1 Mpc. This much larger volume would encompass 
~ 3000 voids with r ^ 5/i _1 Mpc, and would detect voids 
with an efficiency better than 30 per cent, providing a sam¬ 
ple of around ~ 1000 voids. With an extended program on 
PFS of 100 nights, it is possible to construct a tomographic 
map covering ~ 200 deg 2 with the same redshift coverage 
and resolution, providing a tenfold increase in volume, and 
therefore, the number of voids (~ 10 4 ). 

For comparison, similar volumes have been explored 
to find voids in low-redshift galaxy positions, although for 
somewhat larger voids. Pan et al. (2012) searched for r > 
10/i _1 Mpc voids in the Sloan Digital Sky Survey Data Re¬ 
lease 7 main galaxy sample (out to z = 0.1), corresponding 
to a volume V « 10 r /i _1 Mpc 3 , finding ~ 1000 voids. Sut¬ 
ter et al. (2012) also found a similar number of voids in the 
SDSS DR7 main galaxy sample (out to z = 0.2) and the 
luminous red galaxy sample (out to z = 0.44). In total, the 
galaxy samples were split into 6 samples covering volumes 
from 10 6 to almost 10 9 /i _1 Mpc (see their table 2). How¬ 
ever, the larger volume samples were covered by brighter, 
more massive galaxies, with larger separations. By z = 0.1, 
the average galaxy separation in DR7 is already larger than 
5/i _1 Mpc, making it difficult to find statistically significant 


small galaxy voids. Sutter et al. (2014a) provided an update 
to this analysis using the Baryon Oscillation Spectroscopic 
Survey Data Release 9 CMASS sample, split into 6 sam¬ 
ples, each covering about 5 x 10 8 h ~ 3 Mpc 3 . This work found 
~ 1500 voids with large radii (> 20/i _1 Mpc). 

5.2 High-Redshift Void Cosmology 

Much of the recent discussion of voids as cosmological probes 
has focused on them as a means to measuring the Alcock- 
Paczynski (AP) parameter, 

Fap = Ft IdaH, (3) 

c 

where Da is the angular diameter distance and the Hubble 
parameter, H : encodes distortions in the line-of-sight (LOS) 
direction. Note that this measurement measures H(z) di¬ 
rectly, rather than an integral as measured by e.g. Type la 
supernovae. In Fig. 8, we show how specific variations in 
cosmological parameters affect the AP parameter, giving a 
qualitative idea of how accurate these measurements must 
be. Specifically, we show changes in the AP parameter with 
5 per cent variations in Q m , 1 or 2 per cent variations in Hk 
per cent, and 10 per cent variations in wq. 

Lavaux & Wandelt (2012) have proposed measuring the 
AP parameter using the anisotropy of stacked voids, in the 
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Figure 8. Fap as a function of redshift for different cosmological 
models, divided by the prediction from our fiducial model. We 
show 5 per cent variations in f2 m , 1 or 2 per cent variations in 
fik per cent, and 10 per cent variations in wq. For all models, the 
value of the Hubble parameter h has been adjusted to keep the 
angular scale of the CMB oscillations 6\ s = d\ s /r s fixed. State- 
of-the-art measurements of Fap have 5 per cent uncertainties, 
measured at z = 0.5 up to z = 2.4 (Blake et al. 2011; Aubourg 
et al. 2014; Beutler et al. 2014; Samushia et al. 2014). 


context of voids identified in galaxy surveys. Even though 
voids should be spherically symmetric when averaged, each 
void will have a certain random asymmetry that will add 
noise to the global stack. Lavaux & Wandelt (2012) showed 
that the uncertainty due to this intrinsic scatter when aver¬ 
aging N voids can be approximated by (their eq. 35): 


CTAP ^ 1 

Tap ~ VN ’ 


( 4 ) 


In terms of cosmological parameters, Fig. 8 shows that it 
requires substantial changes in cosmological parameter val¬ 
ues, by today’s standards, to produce one or two per cent 
changes in Fap- Therefore, in order to be competitive with 
other cosmological probes, the stack should be done using 
several thousand voids. 

In order to accurately estimate the uncertainty when 
stacking N voids identified in the flux field, we would have 
to study the intrinsic scatter in the asymmetry of the void 
flux profiles, as a function of redshift and void size. We would 
also have to take into account the effect of potential system- 
atics like errors in centering and measuring radii, as well as 
different sources of contamination in the Lya flux. But as¬ 
suming that we would also need ten thousand voids to have a 
one per cent measurement, we can use the discussion above 
to estimate that we could achieve this uncertainty with a 
CLAMATO-like survey over 100 square degrees, or with a 
PFS-like survey over 200 square degrees. 

Moreover, Fap can also be robustly constrained from 
anisotropic measurements of the Baryon Acoustic Oscilla¬ 
tion (BAO) scale (Eisenstein & Hu 1998; Seo & Eisenstein 
2003). BAO measurements typically report ratios of separa¬ 
tions with respect to a fiducial model along the line of sight 
(ay =L cry) and transverse (aj_ =L a±) directions, as well as 
their correlation coefficient (r). One can translate these val¬ 


ues into a ratio of Fap with respect to Fap in the fiducial 
model: 


r _ Fap _ ot± 

with an uncertainty given by: 

g / _ °f . A o rg lP_L 

/ 9 o i o " 

ap a| a a\\a± 


( 5 ) 

(6) 


For instance, recent BAO measurements from the BOSS col¬ 
laboration (Aubourg et al. 2014; Samushia et al. 2014) can 
be translated into ~ 5 per cent measurements of /ap both at 
z = 0.57 (from the galaxy survey) and at z = 2.4 (from the 
Lya survey), raising the bar for measurements from voids. 


6 CONCLUSIONS 

In this paper, we characterized the signal of cosmological 
voids in the high-redshift matter density field and demon¬ 
strated how we can use Lya forest tomographic maps 
to find high-redshift voids. We used a simple spherical 
over/underdensity approach to identifying voids in a large 
cosmological simulation (with a box size of 256 h~ x Mpc or a 
volume of 1.7 x 10 7 h -3 Mpc 3 ) at z = 2.5, resulting in a cat¬ 
alog of - 16, 000 voids with radii of 2 - 12 h -1 Mpc. We also 
tested finding voids with a watershed approach and found 
that the resulting catalog was similar to that produced by 
the spherical overdensity method, but with more complex 
geometries that changed the void centroid non-trivially. For 
simplicity, we used the spherical overdensity void finding 
method throughout. This makes our results somewhat con¬ 
servative, i.e. it is likely that more sophisticated void-finding 
methods will have improved performance. 

Overall, the signature of high-redshift voids in flux is 
similar to what has been found for low-redshift voids in den¬ 
sity. The radial density profile of voids is low (p/p = 0.2 - 
0.4) and rises more steeply closer to the radius of the void. 
One difference we noticed is that the high-redshift voids are 
typically less evacuated than their low-redshift analogues, 
giving them a steeper inner profile and less pronounced 
rise at the edge. The shape of the density profile is clearly 
mirrored in flux with high transmission inside the radius 
(5f — 0.25 - 0.15), and dropping down to the mean flux 
beyond the radius. Interestingly, the radial velocity profiles 
show very little scatter and the mean radial velocity profile 
matches up to the linear theory prediction very well. This 
could be a promising testbed for any (modified gravity) the¬ 
ory predicting differences in void outflow velocities. 

Using our void finding method, we identified voids in an 
ideal flux field and in three tomographic flux maps gener¬ 
ated from mock surveys with spatial samplings of d± = 2.5, 
4, and 6/i _1 (hires, midres, and lores maps). We compared 
the flux void catalogs to the density void catalogs by consid¬ 
ering how well ‘matched’ pairs of voids are in terms of their 
centers and radii. We found excellent agreement between the 
density and ideal flux void catalogs, where 99 per cent of the 
large voids (r > 6/i _1 Mpc) are well-matched. The noise in 
the tomographic maps clearly impacts the efficiency of find¬ 
ing voids, reducing the fraction of well-matched large voids 
down to 66, 48, and 27 per cent in the hires, midres, and 
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lores maps, respectively. However, when we inspected indi¬ 
vidual cases of poorly matched voids, we found that many 
of these are due to noise in the maps artificially breaking up 
or merging high-transmission regions. It is clear that a more 
sophisticated void finder, especially one that models a noise 
component, would perform much better on the tomographic 
maps. Implementing such a method is beyond the scope of 
the current work. 

Using these matching results, we can provide a conser¬ 
vative forecast for the number of voids that can be found 
in dense Lya surveys. Our hires map has a signal-to-noise 
ratio distribution and sightline spacing similar to the ongo¬ 
ing CLAMATO survey. With a sky coverage of one deg 2 , 
the CLAMATO data would produce a tomographic map 
covering V ~ 10 6 /i _1 Mpc 3 , and our proposed void-finding 
method would identify about 100 voids with r > 5/i -1 Mpc 
in such a volume. With a 16 deg 2 survey on the PFS, we 
would identify about 1000 voids with r > 5/i _1 Mpc, al¬ 
though at a degraded purity. A 100-night dedicated Lya 
forest survey across 200 deg 2 on the PFS would increase this 
number by a further order of magnitude to > 10 4 voids. 

These populations of high-redshift voids could be use¬ 
ful for many purposes, including tests of modified gravity, 
as an AP test and for studying high-redshift void galax¬ 
ies. Previous works have considered voids as a clean envi¬ 
ronment for studying galaxy evolution, where galaxies are 
very isolated and their evolution is not complicated by envi¬ 
ronmental effects (e.g. see van de Weygaert & Platen 2011, 
sec. 5). However, existing studies of void galaxies are concen¬ 
trated at low redshift, where such objects are much easier to 
find (van de Weygaert et al. 2011). At low redshifts, the evi¬ 
dence points to the different properties of void galaxies being 
caused by their low stellar mass, independent of other influ¬ 
ence from their void environment (Hoyle et al. 2005; Kreckel, 
Joung & Cen 2011; Tinker et al. 2008). It would be very in¬ 
teresting to see whether similar behavior is seen at higher 
redshifts, where we expect the processes of galaxy forma¬ 
tion could be different. Current galaxy redshift surveys can 
probe only down to L ~ L* in galaxy luminosity at these 
redshifts, and we would naively expect high-redshift voids 
identified through Ly<a forest tomography to also be void of 
such bright galaxies. However, the James Webb Space Tele¬ 
scope and its NIRSPEC spectrograph 3 will have the ability 
to target L ~ 0.3L* galaxies within voids identified through 
CLAMATO and PFS. 

With dense Lya forest surveys covering larger volumes, 
such as a dedicated program on the PFS covering 200 deg 2 , it 
is possible to identify a population of 10 4 voids. Such a large 
number of voids would naively translate to a one per cent AP 
measurement, although this is just a statistical estimate and 
it is possible that there would be larger systematic errors in 
such a measurement. 

As described in App. B, we are providing public access 
to the data used in this project including gridded simulation 
quantities, the tomographic flux maps, a grid of hires flux 
skewers, FoF halo catalogs, void catalogs, and the protoclus¬ 
ter catalog from (Stark et al. 2014). 

We thank Kathryn Kreckel for useful discussions. The 
simulation, mock surveys, and reconstructions discussed in 


3 http://www.stsci.edu/jwst/instruments/nirspec 


this work were performed on the Edison Cray XC30 sys¬ 
tem at the National Energy Research Scientific Computing 
Center, a DOE Office of Science User Facility supported by 
the Office of Science of the U.S. Department of Energy un¬ 
der Contract No. DE-AC02-05CH11231. This research has 
made use of NASA’s Astrophysics Data System and of the 
astro-ph preprint archive at arXiv.org. 
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Figure Al. A density slice centered on a large void with a 
radius of 9.6h -1 Mpc. Just as in Fig. 1, the slice is 40/i _1 Mpc 
across and 6h -1 Mpc into the page. The panels on the left show 
the view along the simulation y axis, while the panels on the 
right show the view along the simultion x axis. (Top) The matter 
density. (Bottom) The corresponding voids from the SO catalog 
(black), the watershed catalog (blue), and the smoothed density 
watershed catalog (red). The dots and dashed line circles show 
the center of the voids and the SO radius or watershed effective 
radius. The triangles show the minimum value (core) points of 
the watersheds. The blue and red colorscales show the projection 
of the points in the watershed (the darker the color, the more 
points into the page). See the text for more details. 


APPENDIX A: COMPARISON OF SO AND 
WATERSHED VOIDS 

In Subsection 2.2, we briefly compared our SO void cata¬ 
log to a set of voids found via a watershed method with the 
same threshold setting. We concluded that it would be more 
straightforward to use the voids found via the SO method, 
mainly due to the complex shapes of the watershed. In this 
Appendix, we show an example void found in both catalogs 
to illustrate this point. In Fig. Al, we show a slice centered 
on a large void. The left and right panels show the same 
void structure from two angles (the xz-plane and yz-plane). 
The top panels show the density field in this region, while 
the bottom panels show the void shape in the different cat¬ 
alogs. The black dot and circle are the center and radius 
of the void found with the SO method. The blue dot and 
circle show the value-weighted centroid and effective radius 
of the void found with the watershed method, and the blue 
triangle shows the ‘core’ point (the minimum value point 
within the shed). We also show the points in the void water¬ 
shed in the bottom panels with the blue colorscale, where 
the color scales with the number of points in the projection. 
In order to damp out some of the complex structure of the 
original watershed void, we also tried running the watershed 
finder on a 2/i _1 Mpc smoothed density field (with an ad- 
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justed threshold of 0.45 p). The watershed void found in the 
smoothed field is shown in red. 

The SO and watershed voids have reasonably similar 
shapes as seen in the yz-plane (right panel). The extent 
of the watershed points (blue region) overlaps almost en¬ 
tirely with the SO circle (black), besides the small wayward 
blue blob at (y,z) = (5,15) h _1 Mpc. However, seen in the 
xz-plane, the voids have very different shapes indeed. The 
slightly overdense region at (x,z) = (—10, 0) h _1 Mpc lim¬ 
its the growth of the SO void, but the watershed region 
reaches around this structure to the underdense region on 
the other side. This extension from the main underdense re¬ 
gion is also seen in the smoothed version of the watershed 
void. The SO radius is 9.6/i _1 Mpc and the watershed effec¬ 
tive radius is 7.8h _1 Mpc, and 9.1 h~ 1 Mpc in the smoothed 
version. Although the radii are all fairly similar, it’s amazing 
to see just how different the extents and centers differ. The 
SO center is 7.2/i _1 Mpc away from the watershed centroid 
and 5.3 h~ 1 Mpc away from the watershed core. At the same 
time, the watershed core and centroid are separated by a 
whopping 11.9/i -1 Mpc, far more than the effective radius. 
Overall, the watershed void spans 39, 26, and 32/i _1 Mpc 
in the x, y, and z directions respectively, meaning the small 
finger-like voids extending from the central underdensity are 
very long. 


APPENDIX B: PUBLIC DATA PRODUCTS 

The data used in this project are available at http: 
//tinyurl.com/lya-tomography-sim-data. We hope that 
making this data publicly available will reduce the barrier 
to future work on Lya forest tomography and high-redshift 
voids. The data release includes gridded simulation quanti¬ 
ties, the tomographic flux maps, a grid of hires flux skewers, 
FoF halo catalogs, void catalogs, and the protocluster cat¬ 
alog from (Stark et al. 2014). Due to space limitations, we 
downsampled the gridded quantities from the full 2560 3 grid 
to a 640 3 grid. Although this process erases some small-scale 
structure, the resolution is still more than enough for our 
purposes. The gridded quantities include the z = 2.5 den¬ 
sity, redshift-space density, flux, real-space flux, and peculiar 
velocities and the z = 0 density and peculiar velocities. We 
also include example Python and C++ sources for reading 
the files. 
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